Level densities of nickel isotopes: microscopic theory versus experiment 
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We apply a spin-projection method to calculate microscopically the level densities of a family 
of nickel isotopes ^®~^*Ni using the shell model Monte Carlo approach in the complete pfgg/2 
shell. Accurate ground-state energies of the odd-mass nickel isotopes, required for the determination 
of excitation energies, are determined using the Green's function method recently introduced to 
circumvent the odd particle-number sign problem. Our results are in excellent agreement with recent 
measurements based on proton evaporation spectra and with level counting data at low excitation 
energies. We also compare our results with neutron resonance data, assuming equilibration of 
parity and a spin-cutoff model for the spin distribution at the neutron binding energy, and find 
good agreement with the exception of ''^Ni. 
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Introduction. Nuclear level densities play an important 
role in determining the elemental abundances in stellar 
nucleosynthesis [Hl^]- However, the microscopic calcula- 
tion of level densities from underlying effective interac- 
tions has been a major theoretical challenge Q. 

The configuration-interaction (CI) shell model ap- 
proach offers an attractive framework for reliable cal- 
culation of nuclear level densities. In this approach, 
both shell effects and correlations are included a priori] 
thereby eliminating the need for a posteriori parameter 
fitting as is done, for example, in empirical approaches 
based on the back-shifted Bethe formula (BBF) 

However, the conventional shell model approach, which 
is based on direct diagonalization of the Hamiltonian in 
a truncated CI shell model space, becomes prohibitively 
difficult in medium-mass and heavy nuclei because of the 
combinatorial increase of the dimensionality of the many- 
body Hilbert space. An alternative is provided by the 
shell model Monte Carlo (SMMC) approach The 
SMMC can be used to calculate thermodynamic observ- 
ables within the CI shell model framework. In partic- 
ular, it has been shown to provide accurate estimates 
of nuclear state densities |13l4l9| . The computational re- 
sources required for SMMC calculations scale gently with 
the size of the single-particle space, enabling calculations 
in very large many-particle model spaces. 

A variety of methods are used to extract level densi- 
ties from experimental data: direct level counting at very 
low excitation energies, neutron or proton resonance data 
at the neutron or proton binding energy, and Ericson's 
fiuctuation analysis [20| at higher excitation energies 
{Ex ^15 MeV in medium-mass nuclei). Level densities 
at intermediate energies can be extracted from charged 



22| 



particle reactions [21[, and using the Oslo method 
in nuclei for which both level counting and neutron res- 
onance data are available. More recently, neutron and 
proton evaporation spectra have been used to determine 
level densities at intermediate energies independently of 
the neutron resonance data [23j. In particular, the level 
densities of a family of nickel isotopes (^^~^''Ni) were ex- 
tracted from proton evaporation spectra [1^. 



The usual density calculated in the SMMC method is 
the state density, in which the spin degeneracy factor 
2 J + 1 of each level with spin J is taken into account. 
However, experiments in the intermediate excitation en- 
ergy regime often measure the level density, in which 
each level is counted only once irrespective of its spin 
degeneracy. In this work we calculate the level densi- 
ties of nickel isotopes (^^-^'^Ni) directly with the SMMC 
method using the spin- proj ection technique described in 
the recent work of Ref. [2^ This method was applied to 
even-mass nuclei only and here we present its first ap- 
plication to odd-mass nuclei. The ground-state energies 
of the odd-mass nickel isotopes are calculated accurately 
using a method recently introduced in Ref. 26. 

We find that our theoretical level densities agree very 
well with the experimental level densities extracted re- 
cently from proton evaporation spectra [23 | and with 
level counting at low excitation energies. We also com- 
pare our SMMC results with level densities extracted 
from neutron resonance data at the neutron binding en- 
ergy, assuming equal densities of positive and negative 
parity levels and a spin cutoff model with rigid-body mo- 
ment of inertia. We find good agreement except for ^'^Ni. 

The shell model Monte Carlo approach. The SMMC 
method is based on a representation of the Gibbs oper- 
ator e^^^ (where H is the CI shell model Hamiltonian 
and /? = 1 /T is the inverse temperature) as a functional 
integral over one-body propagators Ua- describing nucle- 
ons moving in external time-dependent auxiliary fields 
a — (t(t). This is formally expressed in terms of the 
Hubbard-Stratonovich (HS) transformation f27^ 



-PH _ 



= / D[a]G^U, 



(1) 



where is a Gaussian factor. 

The thermal average of an observable O can be written 

as 



J D[a]W^<i>„0„ 



(2) 



where Wa = GcrlTr Uc,\ is a positive definite weight, 
<I>cr = Tr Ua- /\Ti Ua \ is the Monte Carlo sign function and 
= Tr{dU„)/TTU^. In SMMC, the r.h.s. of Eq. @ 
is calculated using stochastic sampling of the auxiliary- 
field configurations a from the distribution Wa- Ther- 
mal averages at fixed number of protons and neutrons 
are calculated in the canonical ensemble with the help of 
particle-number projection. 

The thermal energy E = E{/3) is calculated from 
E{(3) = (H) and the partition function Z{(3) is obtained 
by integrating the thermodynamic relation —d\nZ/dj3 = 
E. The entropy and heat capacity are calculated from 
S = hYZ + I3E and C = -(3'^dE/d(i, respectively. The 
average state density (in which each level with spin J 
is counted 2 J + 1 times) is then obtained by using the 
saddle-point approximation to the inverse Laplace trans- 

1/2 

form of the partition function [28|, p « (27r/?~^C) e'^. 

Experiments, however, often measure the level density 
p (in which each level is counted only once irrespective 
of its spin degeneracy) rather than the state density. It 
was recently shown in Ref. [2^ that level densities can be 
obtained directly within SMMC by calculating the den- 
sity pm of states with given magnetic quantum number 
M for the minimal absolute value of M , i.e., 
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for even-mass nuclei 
for odd-mass nuclei 



(3) 



The density pM can be calculated by projection on 
the total spin component M [l^. Using the saddle- 
point approximation for the M-projected density, we 
find (in analogy with the total state density) pM ~ 

(27r/3-2CM) ^ e'^", where Sm and Cm are, respec- 
tively, the M-projected entropy and heat capacity. In 
this work we use this method to calculate level densities. 
The M-projected heat capacity Cm is calculated inside 
the HS path integral by the method of Ref. [2^ taking 
into account correlated errors and thus reducing signifi- 
cantly the statistical errors in the numerical derivative of 
Em- 

We use the CI shell model Hamiltonian described in 
Ref. [13I in the complete pfg9/2 shell. This Hamiltonian is 
known to provide a good description of the statistical and 
collective properties of medium-mass nuclei in the iron 
region. Its interaction includes the dominating collective 
components of realistic nuclear effective interactions , 
yet has the advantage that it satisfies the modified sign 
rule 0, and is therefore free from the Monte Carlo sign 
problem for even-even nuclei. 

Accurate determination of ground-state energies- In 
order to compare the SMMC level densities with experi- 
mental data, it is necessary to express them as a function 
of the excitation energy. This requires an accurate deter- 
mination of the ground-state energy Eq- 

For even-even medium-mass nuclei there are two accu- 
rate methods to calculate Eq. The first method is based 
on Eq being the limiting value of E{[3) for /? — ;> cxo. At 
large /3, the matrix representing the propagator Ua- in 
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FIG. 1: The parameter Eo (solid circles) versus /3 for ^"Ni in 
the two-level model (see text). The solid line is our estimate 
for the ground-state energy obtained by averaging Eq over the 
/? values within the length of the solid line. The dashed lines 
describe the statistical error in this estimate. 



the single-particle space becomes ill-defined, and we use 
the stabilization method of Ref. [3l| in the framework of 
the canonical ensemble. We can then calculate E{13) for 
several large but finite values of /3 and take an average. 
The second method is based on the two-level model [ij] . 
At large /3, the thermodynamic observables of the even- 
even nucleus are well- approximated by considering just 
the lowest two levels: the 0+ ground state and the low- 
est excited 2+ level. In this model the thermal energy 

E{l3) ~ {H) and (J )^ (where J is the total nuclear spin) 
are (at a given value of /3) functions of Eq and the exci- 
tation energy E^^ of the lowest 2+ state. We can then 
solve these relations to determine Eq and E^ from the 
calculated SMMC values of £;(^) and (J^)^. The final es- 
timates for Eq and i?^^ are then obtained by taking their 
average over moderate to large values of /? (for which the 
two-level model is valid). 

We demonstrate the second method in Fig. [T] for ^°Ni, 
where the parameter Eq is shown as a function of /? (solid 
circles). The final estimate for the ground-state energy 
is the solid line obtained by averaging Eq over the range 
of /3 values spanned by the solid line. The dashed lines 
describe the statistical error of this estimate. 

The projection on an odd number of particles leads to 
a new sign problem at low temperatures, even for the 
good-sign interaction used in this work. This makes it 
impractical to estimate the ground-state energy of an 
even-odd nucleus from direct SMMC calculations. It was 
shown recently in Ref. [2^ that this sign problem can be 
circumvented for the ground-state energy by exploiting 
the asymptotic properties of the single-particle Green's 
functions of the neighboring even-even nuclei. 

For a rotationally invariant and time-independent 
Hamiltonian, the scalar imaginary-time single-particle 



Green's functions are defined as 



Tr^ 



Tr^ e 



(4) 



where v = {nlj) labels the nucleon single-particle orbital 
with radial quantum number n, orbital angular momen- 
tum I and total spin j. Here T denotes time ordering 
and a^miT) = e'^^a^mS"^^ is an annihilation operator 
of a nucleon at imaginary time t(— /3<r</3)ina 
single-particle state with orbital v and magnetic quan- 
tum number m {—j < m < j). Using the HS transfor- 
mation, the Green's functions in Eq. (|4]) can be written 
in a form suitable for SMMC calculations [26l|. 

In the limit /3 — ?> oo and |t| — ^ oo while |t| <^ /3, the 
Green's function for an even-even nucleus with a spin 
J — ground state is well approximated by a single ex- 
ponential, 
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-AB,|r| 



(5) 
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FIG. 2: Ground-state energy estimate Eg versus /3 for ^^Ni, 
obtained from the Green's functions of ^''Ni (solid circles) and 
of "^Ni (open circles). Our final estimate for the ground-state 
energy is described by the solid and dashed lines as in Fig. [T] 



where AEj is the energy difference between the lowest 
spin ,/ = j level of the relevant odd-mass nucleus and 
the ground state of the neighboring even-even nucleus. 
In this asymptotic regime for r, we can calculate A.Ej 
from the slope of In Gi, (r) . Minimizing AEj over all pos- 
sible values of j, we determine the difference between 
the ground-state energy of the odd-mass nucleus and the 
ground-state energy of the even-even nucleus as well as 
the ground-state spin of the odd-mass nucleus. 

For an even-even nucleus with a good-sign interaction, 
Gjy(T), and hence AEj, can be calculated within SMMC 
free of any sign problem. Then, knowing the ground-state 
energy of the even-even nucleus (which can be calculated 
using the method described above) , the ground-state en- 
ergy of the neighboring odd-mass nucleus can be calcu- 
lated accurately. 

In a shell model Hamiltonian with well separated en- 
ergy levels, only a few transitions between the even- 
even and odd-mass nuclei contribute significantly to the 
Green's functions. In this case, the asymptotic region 
can be accessed even at moderate values of /3. 

We demonstrate the method for ^^Ni in Fig. ^ where 
we show the ground-state energy Eq versus (3 for ^^Ni. 
Values of Eg can be determined from both the Green's 
functions of ^"Ni (solid circles) and of ^^Ni (open circles) 
by using Eq. ([5]). The final estimate for the ground-state 
energy of the odd-mass nucleus is obtained by averaging 
over a range of /3 values and is described by the solid 
and dashed lines (see caption of Fig. [2]). The statistical 
errors in this method are comparable to those of the even- 
even nuclei and are an order of magnitude smaller than 
the statistical errors of odd-mass nuclei in direct SMMC 
calculations. 

Table [J summarizes our results for the ground-state en- 
ergies Eq of both the even -mass and odd- mass nickel iso- 
topes. Also shown for the even-mass isotopes are the ex- 
tracted excitation energies of the lowest 2+ levels. They 







Ex 


Exp. 




(MeV) 


(MeV) 


(MeV) 


59Ni 


-251.73(1) 








-259.87(3) 


1.20(7) 


1.33 




-264.38(1) 






62 Ni 


-271.84(3) 


1.29(8) 


1.17 




-275.67(1) 






64Ni 


-282.39(3) 


1.4(2) 


1.35 



TABLE I: SMMC ground-state energies Eo (second column) 
for the corresponding nickel isotopes shown in the first col- 
umn. The third and fourth columns show the theoretical and 
experimental excitation energies E^ of the first excited 2 
states in the even-mass nickel isotopes. The theoretical values 
of Eo and E^ for the even-mass nickel isotopes were deter- 
mined from the SMMC calculations of E{l3) and {J'^)p using 
the two-level model [ij and taking an average over a suitable 
range of j3 values. 



are in reasonable agreement with the experimental val- 
ues. 

Level densities of nickel isotopes. The level densities 
for 59-64]v[i shown Fig. H The SMMC level densities 
(solid circles with statistical error bars) are calculated 
using Eq. ([3]) . The excitation energies are calculated from 
Ex = E — Eq using the ground-state energies Eq listed 
in Table m 

We compare our theoretical level densities with exper- 
imental level densities extracted from different data sets: 
(i) direct level counting (solid histograms) [7] , (ii) proton 
evaporation spectra (open squares) [24| and (iii) neutron 
resonance data (triangle) [7]. 

Level counting is carried out up to a nucleus-dependent 
cutoff energy below which a complete set of energy levels 
is available from spectroscopic experiments. In Table |TT] 
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FIG. 3: Level densities for ^^~^*Ni isotopes: theory versus experiment. The SMMC level densities (solid circles) are compared 
with level counting data (solid histograms) f?!, and with level densities extracted from proton evaporation spectra (open 
squares) [Ij, and (when available) neutron resonance data (triangle) 



we list this cutoff energy and corresponding total number 
of levels used in the direct level counting for each of the 
nuclei. The level densities were obtained using bins of 
size ~ 0.8 MeV. 





No. of levels 
(complete set) 


Cutoff energy 
(MeV) 


En 

(MeV) 


D 

(eV) 


59Ni 


27 


2.713 


8.999 


13.4(9) 


60Ni 


48 


4.613 


11.388 


2.0(7) 


•^iNi 


36 


2.905 


7.820 


13.8(9) 




49 


4.503 


10.596 


2.1(2) 


63Ni 


23 


2.697 


6.838 


16(3) 


64Ni 


49 


4.762 







TABLE IL Total number of levels used for level counting (sec- 
ond column) in a complete set of experimental levels up to a 
certain cutoff energy (third column) . The neutron resonance 
energy En and the mean s-wave resonance spacing 0] are 
shown in the fourth and fifth columns, respectively. 

Recently the Ohio University group [1^ extracted level 
densities of the above nickel isotopes in the intermediate 
energy regime from the measurements of proton evap- 
oration spectra in ^'^Li induced reactions on 54,56,58pg 
These experimental level densities are normalized using 
the level counting data and are shown by open squares 
in Fig. El 

The level density at the neutron binding energy E^i is 



obtained from the mean spacing D of s-wave resonances 
(the values of En and D, where available, are shown in 
Table . The conservation of spin and parity implies 
that only levels with certain spins and parities contribute 
to D. Therefore, to convert the neutron resonance data 
to a total level density it is necessary to make certain as- 
sumptions regarding the distributions of spin and parity. 
We make the usual assumptions that positive and neg- 
ative parity levels contribute equally, and that the s pin 
distribution is described by the spin-cutoff model [29| 
with rigid-body moment of inertia. The corresponding 
level densities at the neutron binding energy are shown 
in Fig. 13] by triangles. 

Overall, our SMMC results are in excellent agreement 
with the experimental level densities over the complete 
experimental energy range. The SMMC level densities 
slightly underestimate the experimental densities in '^^Ni 
and ^''Ni. In ^"^Ni, our calculations are in close agreement 
with the level density extracted from the proton evapora- 
tion data, but are below the level density extracted from 
the neutron resonance data. We note, however, that the 
latter also differs from the proton evaporation results, as 
was already discussed in Ref. [IJ. The reason for this dis- 
crepancy might be that the assumption of parity equili- 
bration, used to extract level densities from the neutron 
resonance data, is not entirely justified at the neutron 
binding energy 

Conclusion. In conclusion, we have presented accu- 



5 



rate microscopic calculations of the total level densities of 
nickel isotopes ^^"^^Ni using the SMMC approach. The 
calculations for the odd-mass nickel isotopes were made 
possible by using two recently developed techniques in 
SMMC. The ground-state energies of the odd-mass iso- 
topes were calculated using the Green's function method 
of Ref. '26!, while the level densities were calculated di- 
rectly using the spin-projection technique described in 
Ref. M 

Our results are in close agreement with experimental 
level densities obtained from level counting and neutron 
resonance data, as well as the more recent level densi- 
ties extracted from proton evaporation spectra. It will 
be interesting to apply the formalism developed in this 
work to other mass and energy regimes that are relevant 



in stellar nucleosynthesis and are not yet accessible in 
experiments. 
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